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Abstract 

We calculate the probability for dilepton production in central relativistic 
heavy ion collisions due to the 77 mechanism. This is a potential background 
to more interesting mechanisms. We find that this mechanism is negligible in 
the CERES experiments. Generally, the contribution due to this mechanism 
is small in the central region, while it can be large for small invariant masses 
and forward or backward rapidities. A simple formula based on the equivalent 
photon approximation and applications to a possible scenario at RHIC are also 
given. 
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I. INTRODUCTION 



In peripheral relativistic heavy ion colhsions there are huge effects due to the strong 
electromagnetic fields, while strong interactions among the colliding ions can be virtually 
neglected, due to their short range In central collisions, on the other hand, the strong 
interactions between the colliding ions completely overwhelm the effects due to the elec- 
tromagnetic interaction. It is the purpose of this paper to study the effects due to the 
interactions of the very strong (coherent) electromagnetic fields in the central collisions, see 
Fig. (|lD. The modification of the fields during the collision due to the strong interactions 
will be neglected for simplicity. We are mainly interested in the order of magnitude of the 
effects. Not unexpectedly, the production of strongly interacting particles (like tcq) via the 
77-mechanism is negligible. On the other hand, the photon-photon mechanism could play 
a role in the production of e^e~-pairs (and also muon pairs). In recent experiments low 
mass electron pairs are measured in p-Be, p-Au, and S-Au collisions. Even heavier systems 
like Au-Pb have recently been studied For the proton induced interactions, the spectra 
are explained essentially by electron pairs from hadronic decays, whereas for the heavier sys- 
tems an enhancement is obtained, which is suggested to come from two-pion annihilation. 
In this paper we are going to calculate the contribution of the 77-mechanism to the electron 
pairs. 

In section ^ the general theoretical framework is presented. Calculations within the 
framework of lowest order QED |^ are given in section JT| and compared to order of mag- 
nitude estimates using the equivalent photon approximation. In these calculations we have 
used the same kinematical constraints as are used in the CERES experiment. We also present 
calculations for a possible scenario at RHIC. Our conclusions are given in section 



II. COLLISIONS OF QUASIREAL PHOTONS IN CENTRAL HEAVY ION 

INTERACTIONS 

A. Lowest order QED approach 

The vector potential of the combined electromagnetic field of two heavy ions, which 
travel with constant four velocities u^^'"^^ and with an impact parameter b between them is 
given by 



+Z2F2iq^)u^;^^5{qu^^^)exp{-tqb/2) 



(1) 



where the four-velocities are given by m(^) = (7, 0, 0, 7/3) and m^^) = (^^ q, 0, -7/?) 0. In the 
following we restrict ourselves to central collisions with b ^ 0. The nuclear charges and form 
factors are denoted by Zi and Fi{q^), respectively {i = 1,2). It should be noted that due 
to the strong interaction between the heavy ions, the nuclear charge distributions (which 
give rise to the electromagnetic field) will change; both stopping, as well as, a spreading out 
will take place. Such effects were considered in . In addition there is also pair production 
from the stopping itself H|. This pair production due to the bremsstrahlung will strongly 
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depend on the actual stopping. We will neglect these stopping effects here. The collision 
of the two electromagnetic fields (equivalent photons) can then give rise to all kind of final 
states /, like e+e~ pairs, ttq, etc.. 

The matrixelement for e^e~ pair production is given by 

M = -ie^u{p_) 

where u{p-) and f are the Dirac spinors describing the produced electron and positron, 
respectively. For impact parameter b equal to zero, this can be calculated completely ana- 
lytically as was done in The integration over all final states was then done with a MC 
integration P,p!0[] in the program code BORNZERO. The form factor of a pointlike charge 
distribution {F{q) = 1) contains unrealistically high Fourier components, which can lead 
to spurious effects. Already the total probability diverges logarithmically with 6 — in 
this case, see 0. In order to get realistic results a nuclear form factor is needed in the 
calculations. We use here either a monopole form factor 
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P'monopole(^Q ) "^^2 ^2' 

where A is chosen to reproduce the correct rms radius of the nucleus, as well as, one which 
is the sum of two monopole form factors. We use them as a check on the sensitivity of our 
results on the detailed form of the form factor. Compared to a "realistic" form factor, the 
monopole form factor does not fall off rapidly enough for large q. On the other hand this 
means, that our result is an upper bound. 



B. Equivalent Photon approximation (EPA) 



The impact parameter dependent equivalent photon number N{uj,p) for a monopole 
form factor is given as |1TT|,|5| 



N{u;,p) 
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A simple approximation neglecting the contribution from the nuclear interior is 
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R < p < 'f/uj 



(5) 



and otherwise. The radius of the ion is given by i? = 1.2 fm A^/\ With this expression a 
simple formula for the product of the equivalent photon numbers integrated over p can be 
given: 
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d^pN{u^,p)N{uj2,p) 



1 - a^e^l^l 



(6) 



where Y = \j1 ln(ti;i/co'2) and a = i?>M/27 < 1, with i?> the larger one of the two nuclear 
radii and M the invariant mass of the produced system. 

In the EPA the transverse momentum of electron and positron are the same but in 
opposite directions (pt(e"'") = —pt{e~)). We have the following expression for the probability 
to emit an e^e~ pair with rapidities ?/+ and y- and transverse momentum pt(e^) 



d^Ptdy+dy^ 



[ d^pN{uJi,p)N{uj2,p)-^ 
J n dt 



(7) 



where da/dt is the cross section for the subprocess 7 + 7 ^ e+ + e~, see, e.g., . 

Similarly the probability to produce a final state with invariant mass M and rapidity Y 
is given by 



d'P ZlZlc? 



dMdY 



(8) 



where cttt is the corresponding 7 + 7 ^ / photon-photon fusion cross section. 
From this equation it can be easily seen that, e.g., ttq production is very small: 



dP 

dy 
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(9) 



For y = we find (for 7 large enough so that a <^1) that the pion production probability 
dP/dY f» 1.3 X 10"^ for Pb-Au or Au-Au colhsions. 

Since the vr" is the lightest hadron, it is also the hadronic state which will be produced 
with the highest probability. Therefore the production of hadronic final states by the elec- 
tromagnetic fields will be negligible compared to the hadronic production. We can therefore 
concentrate in the following only on e^e~ pairs, similar calculations are of course also pos- 
sible for muon or r pairs. As there was some concern regarding e^e~ pairs coming from the 
electromagnetic production in connection with the experiments 
exactly. 
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III. NUMERICAL RESULTS AND COMPARISON WITH EXPERIMENT 

We want to compare our results directly with the measurements of the CERES exper- 
iment and also want to make a prediction for a possible scenario at RHIC |1^. For 
comparison with CERES we have used the same kinematical restriction {pt > 200MeV for 
both S-Au and Pb-Au collisions). The CERES experiment also has a restriction on the 
relative angles between electron and positron (Bee > 35mrad). This was not taken into 
account, as tests have shown that — apart from the very small invariant masses — this 
condition was always fulfilled. 

The results of our calculations are shown in Figs. We compare directly 

{d'^Nee/ dydruee) / {dNch/ dy) , using dNch/dy = 125 for S-Au colhsions, dN^h/dy = 235 (the 
average of two experimental runs of CERES for Pb-Au collisions.) We show the quantity 
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(PN/dydrriee, which is defined to be the differential probabihty with respect to the y of 
either the electron or the positron and integrated over the allowed y range of the other 
{2.1 < y < 2.65 in the lab frame, corresponding to —0.93 < y < —0.38 in the "center of 
velocity" frame), following 0. 



For the RHIC scenario, see Fig. ^, we have used the conditions as discussed in []T^ 
Pt > 200 MeV, -0.35 <y < 0.35. dNch/dy was assumed to be 1100. 

We note the following points: The results of our calculations are well below the experi- 
mental results of and 0]. Since the cross-section for the 77 mechanism scales with ZfZ2 
we expect the strongest effects for the heaviest systems. For the p-Be and p-Au case at 
450 GeV 0, the 77 mechanism is only a small fraction, as expected. We do not need to 
show this. Even for the heaviest systems 0], the effects are well below the measurements. 
Therefore we can safely conclude that the 77-mechanism plays only an insignificant role for 
the experimental results of 

Calculations with both the single and the double monopole form factor are shown for Pb- 
Au and Au-Au collisions. The results for the CERES experiment are found to be sensitive 
to the form of the factor both at low and high invariant masses, those for RHIC only for 
the low- invariant mass region {rriee < 2pt). The sensitivity in this range is rather easy to 
understand: In order to be in a low invariant mass state, even though the momenta is at 
least Pt, electron and positron cannot be produced in opposite transverse directions, as is 
assumed in the equivalent photon approximation. In the EPA the invariant mass of the 
lepton pair is given by 

mee = 2px cosh [1/2 {y+ - . (10) 

Therefore the sum of their momenta is larger than the typical range of the form factor (about 
50 MeV for Pb or Au). This makes them sensitive to the high q part of the form factor. 

The discrepancy at higher invariant masses at the CERES experiment comes from a 
different reasoning: In order to produce a pair with, for example, rriee = 1 GeV and with 
Y = 0, one needs photons with energy about u = 500MeV^ (estimated from the equivalent 
photon approximation). The kinematics of this process relates this with the virtuality of 
the photon: 

Q' = ql + {uhf . (11) 

With 7 ?^ 10 the virtuality of the photon is then at least 50 MeV at CERES, but only 
at least 7 MeV at RHIC. Therefore at CERES one is again sensitive to the form factor at 
large virtuality. 

For an order of magnitude estimate, we can use the EPA. We use Eq. (|^ together with 
the equivalent photon number of Eq. (|^) and the exact expression for da/dt, see, e.g., 
Tsj , and integrate over d'^pt. 

One can find also a simple analytic formula for the probability to produce dileptons in 
the central region. We can neglect the square of the electron mass as compared to the values 
of the Mandelstam variables s and t of the subprocess 7 + 7 — > e"*" + e~ and find 



da lixa^ (s + t t 
di P \ i s + t 



. > . +^)- (12) 
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For y_ we have i = 



— s/2 and we find 



In Figure ^ the different approaches are compared both for Au-Au colhsions at RHIC 
(-0.35 <y < 0.35, Ay = 0.7) and Pb-Au colhsions at CERES {2.1 < y < 2.65, Ay = 0.55). 
One finds that EPA agrees quite weU at invariant masses above 2pt with the exact calcula- 
tion. The analytic expression given above is in good agreement with the results at RHIC, 
but as explained above, is not valid for the conditions at CERES, as the parameter a is larger 



than 1 (EPA would predict in this case (a > 1) a value of zero!). Equation ([T2|) can be rewrit- 
ten in terms of a differential cross section da/dQ, which results in a [1 + cos^(^)] / sin^(6') 
dependence on the scattering angle in the center of mass system. From this one sees that 
the importance of the 77 mechanism, when one goes to more forward or backward angles 
(larger rapidities), increases as compared to the other mechanisms for dilepton production, 
see also M. 



IV. CONCLUSIONS 

A comparison with the measured e~^e~ pairs in central collisions at CERES shows, that 
pairs produced electromagnetically from the strong electric fields, are only a small back- 
ground. We think that this calculation has, for the first time, convincingly shown that the 
77 pairs are negligible. This is mainly due to the experimental conditions on large trans- 
verse momenta and large invariant masses. The majority of the electron-positron pairs is 
produced with a invariant mass of the order of 1 — 10 rrie with rather small pt. Also the 
photon-photon mechanism is strongly peaked at very forward angles and therefore 77-pairs 
will show up prominently under such kinematical conditions 

We extended our calculations also to RHIC energies. In this energy regime, 77-pairs 
are produced with invariant masses up to several GeV. Therefore the 77 mechanism can be 
potentially dangerous in this case. Again, our calculations show that this is only a small 
effect in the central region, where we compare our calculations with the predictions given in 
1^. In any case, with the present calculations the contributions due to the 77 mechanism 



are well under control. They are reassuringly small in the central region, where one is looking 
for the effects of hot and dense matter. 
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FIGURES 




FIG. 1. The strong electromagnetic fields surrounding the heavy ions are also present in central 
collisions (6 = 0) and will lead to pair production processes due to the 77 fusion mechanism. 
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FIG. 2. d?Nee/dydmee is shown for the electromagnetic electron-positron pair production for 
central S-Au collisions at CERES {pt,min > 200MeV/c); also shown are the experimental results of 
the CERES experiment O]. 
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FIG. 3. (fNee/dydmee is shown for the electromagnetic electron-positron pair production for 
central Pb-Au collisions at CERES {pt^min > 200MeV/c) and for two different form factors; see 
text for details. 
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FIG. 4. Same as Fig. ^ for central Au-Au collisions at RHIC; see text for details. 
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FIG. 5. Different approximations to the pair production probability are compared. Shown are 
the results of the full calculation (BornZero), the full EPA calculation (EPA, Eq. (^) with of 
Eq. (^) and the analytic expression of Eq.([T^) (analytic). 
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